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Abstract 

We study low-energy nucleon Compton scattering in the framework of baryon chiral perturbation 
theory (B%PT) with pion, nucleon, and A (1232) degrees of freedom, up to and including the next- 
to-next-to- leading order (NNLO). We include the effects of order p 2 , p 3 and p 4 /A, with A ~ 300 
MeV the A-resonance excitation energy. These are all "predictive" powers in the sense that no 
unknown low-energy constants enter until at least one order higher (i.e, p 4 ). Estimating the 
theoretical uncertainty on the basis of natural size for p 4 effects, we find that uncertainty of such 
a NNLO result is comparable to the uncertainty of the present experimental data for low-energy 
Compton scattering. We find an excellent agreement with the experimental cross section data up 
to at least the pion-production threshold. Nevertheless, for the proton's magnetic polarizability 
we obtain a value of (4.0 ± 0.7) x 10 -4 fm 3 , in significant disagreement with the current PDG 
value. Unlike the previous ;\PP studies of Compton scattering, we perform the calculations in 
a manifestly Lorentz-covariant fashion, refraining from the heavy-baryon (HB) expansion. The 
difference between the lowest order HB^PT and B%PT results for polarizabilities is found to be 
appreciable. We discuss the chiral behavior of proton polarizabilities in both HBxPT and B%PT 
with the hope to confront it with lattice QCD calculations in a near future. In studying some of 
the polarized observables, we identify the regime where their naive low-energy expansion begins to 
break down, thus addressing the forthcoming precision measurements at the HIGS facility. 

PACS numbers: 13.60.Fz - Elastic and Compton scattering, 14.20.Dh - Protons and neutrons, 25.20.Dc - 
Photon absorption and scattering, 11.55.Hx Sum rules 
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I. INTRODUCTION 

Compton scattering off nucleons has a long and exciting history, see Refs [1-3] for recent 
reviews. The 90's witnessed a breakthrough in experimental techniques which led to a 
series of precision measurements of Compton scattering [4-9] with the aim to determine the 
nucleon polarizabilities [10, 11]. 

Many theoretical approaches have been tried in the description of nucleon polarizabili- 
ties and low-energy Compton scattering. The more prominent examples include dispersion 
relations [12-18], effective-Lagrangian models [19-22], constituent quark model [23], and 
chiral-soliton type of models [24-28]. There has been as well a significant recent progress 
in approaching the subject from first principles — lattice QCD (1QCD). The present 1QCD 
studies are based on the external electromagnetic field method [29, 30], and even though the 
actual results for the nucleon have been obtained only in quenched approximation, the pion 
and kaon polarizabilities have been calculated with dynamical quarks [31]. The full-lQCD 
calculations for the nucleon will hopefully be done in a near future. 

In this work we exploit another theoretical approach rooted in QCD, namely, chiral per- 
turbation theory (%PT) [32-35]. The very first xPT calculation of nucleon polarizabilities, 
published in 1991 by Bernard, Kaiser and Meifiner [36], quotes the result shown in the <^(p 3 ) 
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column of Table I. In the same column, the numbers in brackets show the result of the so- 
called heavy-baryon (HB) expansion [38]. Here it means that one additionally expands the 
full result [36] in powers of m^/M^v, the ratio of the pion and nucleon masses, and drops all 
but the leading terms (cf. Appendix A). The 0(p 3 ) HB;\TT result thus corresponds to the 
static nucleon approximation. The relativistic effects are systematically included in HB%PT 
at higher orders, but nonetheless even leading order HB%PT result is widely considered to 
be more consistent than the B%PT (i.e., fully relativistic) one. The reason for that is that 
the full relativistic evaluation of the chiral loops may yield contributions which are of lower 
order than is given by the power-counting argument. This "pathology", however, does not 
arise in the case of polarizabilities at 0(p 3 ), so as far as power counting is concerned, the 
B%PT result is as good here as the one of HB%PT. But even more generally, chiral symmetry 
ensures that the power-counting violating terms to be always accompanied by low-energy 
constants (LECs), hence they can simply be removed in the course of renormalization of 
those LECs [39, 40]. In simpler terms, there is no problem with power counting in B^PT. 

The present state-of-the-art studies based on pion and nucleon degrees of free- 

dom [41, 42] utilize the HB expansion. They find, however, that despite the very reasonable 
values for polarizabilities, the &{p z ) and even 0(p 4 ) results for the Compton-scattering cross 
sections show significant discrepancy with experimental data starting from energies of about 
120 MeV, especially at backward kinematics. The inclusion of the A(1232)-resonance as an 
explicit degree of freedom helps to remedy this discrepancy in the cross sections [43, 44]. 
However, it comes at an expense of a large contribution to the polarizabilities [45]. This 
A-contribution is highly unwanted in HB%PT, since polarizabilties come out nearly perfect 
already in the theory without the A (cf. the numbers in brackets in Table I). There is 
no natural solution to this problem. One is bound to either omit some of the A contribu- 
tions by "demoting" them to higher orders [43], or cancel them by "promoting" some of the 
low-energy constants (LECs) to lower orders [44]. 

Such an apparent failure of x?T is sometimes attributed to certain "cr-meson" contri- 
butions [46], which xPT misses. Of course, while the a-meson of the linear sigma model 
is included in x?T, the contribution from the /o(600) is not, but it is doubtful that the fo 
can explain it; its two-photon coupling is too small. Alternatively, studies based on disper- 
sion relations suggest that some essentially relativistic effects, discarded in HB^PT as being 
higher order, are in fact important because of the proximity of cuts in both pion mass and 
energy [47-49]. 

In our present study we verify the latter scenario and perform the calculations in a 
manifestly Lorentz-covariant fashion, refraining from the use of the heavy-baryon formalism. 
The HB^PT results can then be recovered by simply expanding in powers of pion mass over 
the baryon mass, itl^/Mb- We thus are coming back to the original (relativistic quantum 
field theory) ways [36]. The difference with the original work [36] is that we compute the 
Compton scattering observables, not only the scalar polarizabilities, and that we include the 
A(1232) in addition to the pion and nucleon degrees of freedom. 



Table I shows the results of both manifest-covariant and HB calculations at all the "pre- 
dictive" orders, i.e., below €?(p 4 ) — the order at which the unknown LECs start to enter. A 
natural estimate of the <ff(p A ) contribution, given in the corresponding column, can serve as 
an error bar on the ;\PT prediction. A detailed discussion of these results can be found in 
Sect. IV. It can be noted, however, how significant the differences are between the exact and 
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BxPT (HBxPT) 


PDG 






^(p3) + ^(p 4 /Z\) 


^(p 4 ) est. 


[37] 




6.8 (12.2) 


10.8 (20.8) 


±0.7 


12.0 ±0.6 


pip) 


-1.8 (1.2) 


4.0 (14.7) 


±0.7 


1.9 ±0.5 



TABLE I: Predictions of baryon xPT for electric (a) and magnetic (/?) polarizabilities of the proton 
in units of 10 -4 fm 3 , compared with the Particle Data Group summary of experimental values. 

the HB results. This is of course not the first and only example where B%PT and HB^PT 
are in dissent, see e.g., the case of 7iV — > A transition [50, 51], or the baryon magnetic 
moments in SU(3) [52, 53]. These differences can often be significantly diminished by slight 
improvements of the HB calculations, such as readjusting the position of the thresholds to 
have them in the exactly correct place [41]. It is not yet clear, however, how to systematically 
derive such improvements from the HB formalism itself. 

As to why the orders considered here are predictive, any chiral power-counting scheme 
will tell us that the expansion of the Compton amplitude begins at order p 2 , and that p 4 
is the order where the first unknown LECs should enter. In between there are p 3 and the 
A-excitation effects. The counting for the latter is itself a subject of controversy related to 
the issue of how to count the A-nucleon mass difference: A = Ma — m 300 MeV. In the 
hierarchy of chiral symmetry breaking scales, A is neither as light as the scale of explicit 
symmetry-breaking, ~ 150 MeV, nor as heavy as the scale of spontaneous symmetry- 
breaking, in f-x ~ 1 GeV. We treat A as an independent light scale with the power-counting 
rules defined in Sect. II. In any case, the leading A effects come before p 4 . 

To recapitulate, in this work we compute the contributions to Compton amplitude up 
to, but not including, ff(p 4 ) in B%PT with A's. This is a complete next-to- next-to-leading 
order (NNLO) calculation which is entirely expressed in terms of only known LECs. The 
details of these calculations are given in Sect. III. Polarizabilities and their chiral behaviors 
are discussed in Sect. IV while the results for observables are shown in Sect. V. 

Some of these results have recently been reported in a letter [54]. The present paper is 
more comprehensive and self-contained. 

II. CHIRAL LAGRANGIANS AND POWER COUNTING 

The method of constructing the chiral SU(2) Lagrangians with pion and nucleon fields 
is well known [34, 35, 55], and the inclusion of the A-isobar fields in a Lorentz-covariant 
fashion has recently been reviewed [56]. We shall list here only the terms relevant to the 
present work. The strong-interaction piece is given by 

f 2 

= J —tr (d^Ud^ + 2B (UM^ + MU^)), (la) 

C$ = N(i$-M N + ^ + g A ^ 5 )N, (lb) 

4 1} = A4tY uX d x -M Al ^)A, + ^\tNT a Y u \d,A u )tT(a x r a ) + R.c.}, (lc) 

where U is the 577(2) pion field in the exponential parameterization: U = exp(in a r a / f), f 
is the pion decay constant in the chiral limit, M is the mass matrix of light quarks, and B 
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is a proportionality factor that can be related with the value of light quark condensate [34] . 
In turn, N denotes the isodoublet Dirac field of the nucleon, is the nucleon mass, and 
qa is the axial-coupling constant, both taken at their chiral-limit value, and the vector and 
axial- vector chiral fields above are defined in terms of the pion field, 7r a (x), as 

v, = \ r a v;{x) = 1 (u + u%u) , (2a) 
a M = \ r a a*(x) = 1 ( u t d^u - u d^) , (2b) 

where u = exp(in a T a /2f) = U 1 ^ 2 . Finally, A u is the Delta isobar Rarita-Schwinger field 
with mass Ma, and Ha is the nNA coupling constant whose value is fixed to the A — > ttN 
decay width of 115 MeV. The antisymmetrized products of Dirac matrices in the above 
equations are defined as: ^ = ^(7^7^ — 7^7^) and 7 MiyA = |(7 Miy 7 A + 7 A 7^). The isospin 
1/2 3/2 transition matrix T is normalized such that T a T bt = \(25 ab - ie abc r c ). 

The electromagnetic interaction is added as usual through the minimal substitution: 

d^N -> d^N - ieA^(l + r 3 )N, (3a) 
d^ a - d,n a - eA^ ab3 7i b , (3b) 

where is the photon field. The minimal coupling of the photon to the Delta field gives 
contributions to Compton scattering which are of higher orders than the ones considered in 
this work. 

There is as well a number of nonminimal terms: 



£ (2) = 

A 2M N (M N + M A ) 



3e 



NT 3 (ig M F^ - g E l 5 F^ d,A u + H.c, (4b) 



e 2 



(4) _ /? F^Vq 

^wzw ~~ 327r 2 / 

Here, and -F M1/ are the photon field strength tensor and its dual tensor defined as 
piiv _ A" — d" A^ , = \e^ vpX F p x] n p (/c n ) stands for the proton's (neutron's) anomalous 
magnetic moment; g& and gM are 7^ A electric and magnetic couplings, respectively, which 
are well known from the analysis of pion-photoproduction P 33 multipoles [51]. Here we differ 
from the strategy adopted in Ref. [43], where, in the absence of any %PT analysis of pion- 
photoproduction at the time, the values of gE and gM were fitted to Compton scattering 
data, with a rather unsatisfactory result. The precise values of all the parameters used in 
the present work are given in Table II. 

Inclusion of the A-isobar fields in a Lorentz-covariant fashion raises the consistency prob- 
lems of higher-spin field theory [57-59]. The A-isobar couplings used here possess the prop- 
erty of invariance under a gauge transformation: A^ — > A^ + <9 M e, where e is an arbitrary 
spinor field. This ensures the decoupling of unphysical spin- 1/2 degrees of freedom and 
eliminates the consistency problems [60, 61]. It is far less straightforward to reconcile this 
extra gauge symmetry with other symmetries of the chiral Lagrangian. For recent progress 
see Refs. [56, 62, 63]. 

Coming to the power counting, for the pion and nucleon contributions we shall use the 
usual scheme [35], i.e., a graph with Vk vertices from & k \ L loops, N n pion and N N nucleon 



5 



G{v 2 ) 


£ = ^ , M N = 938.3 MeV, he = 


197 MeV- fm 




g A = 1.267, U = 92.4 MeV, = 


139 MeV, m w o = 136 MeV, n p = 1.79 




M A = 1232 MeV, h A = 2.85, g M = 


= 2.97, g E = -1.0 


0(p*) 


2 

a ,Po - ± 4nM % 





TABLE II: Parameters (fundamental and low-energy constants) at the order they first appear, 
lines is of order p n with 



n 



J2 kV k + AL- 2N W - N N . 



(5) 



In the case of Compton scattering by p we understand the photon energy or/and the pion 
mass as compared with 47rf n ~ 1 GeV. 

The graphs with As are more tricky because for small p they go as 



S A 



1 



p±A 

rather than simply 1/p as the nucleon propagators. The new scale 

A = M A - M N ~ 293 MeV 



(6) 



(7) 



is neither as light as m % nor as as heavy as 4:iif n , hence can and will be treated independently. 
For energies comparable to the pion mass we choose to additionally expand in p/A, and hence 
the A propagator counts as I /A, while a graph with A^a internal lines contributes to order 



p" 



N A 



(8) 



For defmiteness, when needed, we count A 2 to be of &(p), i.e., the "<5 counting" scheme [43]. 

For the power counting in the region where the energies are of order of A (the resonance 
region) see [43, 51]. Hereby we limit ourselves to the low-energy region, 



p ~ m w < A < 4?!-/;. 



(9) 



III. COMPTON AMPLITUDE AT NNLO 
A. Graphs and the nucleon field redefinition 

The chiral expansion for the Compton amplitude begins with graph (1) in Fig. 1 and its 
crossed counterpart. Nominally they both are of <ff{p)-, however, together they simply give 
the Thomson amplitude, which is of &{p 2 ). We thus refer to ff{p 2 ) as the leading order 
(LO). The other graphs in Fig. 1 contribute to &(p 3 ), the next-to-leading order (NLO). 
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(2) 



(3) 



(i) 



FIG. 1: The tree graphs evaluated in this work. Graphs obtained from these by crossing and time 
reversal are not shown, but are evaluated too. Dots stand for the (leading) jNN and irNN vertices 
from £W, whereas crossed circle denotes the coupling via anomalous magnetic moment. Filled 
square stands for the WZW-anomaly 7r°77 vertex from 





FIG. 2: The one-loop graphs contributing to Compton scattering at 0(p 3 ). Graphs obtained from 
these by crossing and time reversal are not shown. 



At NLO we also have the one-loop contributions shown in Fig. 2, but before evaluating 
them we make a redefinition of the nucleon field, N — > £N, where 



£ = exp 



if) A 7T a T a 



2/ 



"75 



(10) 



The first-order chiral Lagrangian Eq. (lb) then becomes: 1 

C' N {1) = Wt(ip-M N + g A t<tiSN 

= N(0- M N ) N + M N N (1 - £ 2 )N 



(11) 



The two Lagrangians are equivalent, in the sense of equivalence theorem, however, may 
have drastically different forms when expanded in the pion field. For the one-loop contri- 
butions to Compton scattering it is sufficient to expand up to the second order in the pion 
field: 

1 



— r a e abc 7r b d^n c + 



e = 1+ 27 T7r 75 8f 



9A n 2 + {n 3y 



(12a) 
(12b) 
(12c) 



1 In our conventions -fl — 75, hence $ = exp(—igir a T a j5/2f 7T ), ££t = 1. Note also that N — > N£, and 
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(5) 



(6) 



(7) 





(9) 



(12) 




(13) 



FIG. 3: The NLO loop graphs evaluated in this work. Graphs obtained from these by crossing and 
time reversal are not shown, but are evaluated too. 



The original and the redefined Lagrangians take, respectively, the following form: 



C 



N 



m - N (t3-M N + ^r a dn a l5 



T a e abc ir b 3ir c I N + 0(tt 



4/ 2 



(13) 



C 



/(i) 

N 



9A , r -2 
2/2 



(9A - I) 2 _ , 



T a e abc ir b 3ir c I N + Gi'K 



4/2 



(14) 



The major difference between the two forms is that the pseudovector nNN coupling is 
transformed into a pseudoscalar one, while the Weinberg-Tomozawa ttttNN term, which 
resembles a p-meson exchange, gets replaced by an isoscalar term akin to the remains of 
an integrated-out a-meson in the linear a model. The isovector ttttNN term, which is now 
proportional to — l) 2 , does not give any contribution to Compton amplitude at one-loop 
level. 

Also, in the NLO loops, the photon couples only minimally, i.e., to the electric charge of 
the pion and nucleon. Now that the pion couples to the nucleon via pseudoscalar coupling, 
there is no Kroll-Ruderman ('ynNN) term arising, and hence the number of one-loop graphs 
is reduced. The resulting expressions for amplitudes become simpler. 

As a result, the loop graphs shown in Fig. 2 with couplings from the Lagrangian Eq. (13) 
transform to the graphs shown in Fig. 3 with the couplings from Eq. (14). We have also 
checked explicitly that the two sets of one-loop diagrams give identical expressions for the 
Compton amplitude. 

Although the main purpose of the above field redefinition is to simplify the calculation, 
it does give more insight about the chiral dynamics. First of all, it explains how Metz and 
Drechsel [64] , calculating polarizabilities in the linear a model with a heavy a-meson, obtain 
to one loop exactly the same result as B%PT at 0(p 3 ) [36]. Secondly, observing that graphs 
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(12) and (13) vanish in the forward kinematics, we can see that chiral symmetry plays less 
of a role in the forward Compton scattering at order p 3 . 

Going to I A) we encounter the graphs with one A-isobar propagator shown in Fig. 4. 
The nucleon-field redefinition does not affect these contributions at this order. 

Note that the graphs where photons couple minimally to A contain more than one A 
propagator and therefore should be suppressed by extra powers of p/A. However, their 
lower-order contributions are important for electromagnetic gauge invariance and therefore 
for the renormalization program. In particular, the lower-order contributions of chiral loops 
should not affect the result of the low-energy theorem (LET) [67], and this condition is 
automatically satisfied for a subclass of graphs which obeys gauge invariance. The loop 
graphs in Fig. 4 form such a subclass for the case of neutral A. In reality the A comes in 
four charge states (isospin 3/2), and hence a gauge invariant set will in addition have the 
higher-order graphs where photon couples minimally to the A. To make the subclass of 
loop graphs in Fig. 4 gauge invariant without the higher-order graphs, we used the following 
procedure: 

- The one-particle-irreducible (1PI) graphs, Fig. 4(15-18) are computed with the correct 
isospin factors, i.e., summing over all charge states of the A. The isospin factors for 
the one-particle reducible (1PR) graphs (19-22) are chosen such that their ratio to the 
isospin factors of 1PI graphs is the same as in the neutral A case. 

This procedure automatically ensures exact gauge invariance and thus effectively includes 
the lower-order contributions of the one-loop graphs with minimal coupling of photons to 
the A. In case when the latter graphs are included explicitly, the isospin factors of 1PR 
graphs can be restored to actual values. This, however, will not affect the result at the order 
considered here. 

The graphs in Fig. 3 and in Fig. 4 were eventually computed by us with the help of the 
symbolic manipulation tool FORM [65] and the LoopTools library [66] using dimensional 
regularization. 




(17) (18) (19) 




FIG. 4: The NNLO graphs evaluated in this work. Double lines denote the propagator of the A. 
Graphs obtained from these by crossing and time reversal are not shown, but are evaluated too. 
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B. Renormalization 



In accordance with the LET [67], the loop contributions shown in Fig. 3 and in Fig. 4 may 
contribute to the renormalization of nucleon mass, field, charge, and anomalous magnetic 
moment. We have adopted the on-mass-shell renormalization scheme, not the extended on- 
mass-shell renormalization (EOMS) [40]. The difference is that in EOMS the above-listed 
quantities are taken to be at their chiral limit value while here we simply take the values at 
the actual pion mass, see, e.g., Table II. 

Let us discuss first the contributions to nucleon self-energy and ^NN vertex correspond- 
ing to the nucleon loops (5-7) in Fig. 3. The corresponding amputated diagrams give 
contributions to nucleon self-energy (5) and to ^NN vertex (6), (7). More specifically, the 
contributions to the self-energy and the '-/NN vertex can be written in the following form: 

iZ(p) = iE(M N )+i(p- M N )E'(M N ) 

+S(p-M N ), (15) 
ir»(p,p s ) = iYF 1 (p 2 s )-iY l/ qvF 2 (p 2 s )+i^F 3 (p 2 s ) 

+i( Ps -M N )^F 4 (p 2 s ), (16) 

where p is the initial nucleon momentum, q is the initial photon momentum, and p s = p + q. 
The function S is finite and its expansion in powers of f — M N starts from a quadratic term. 
Of all the functions F 1 . . . F 4 that contribute to the >yNN vertex only Fi is divergent. It 
contributes to the renormalization of charge. Function F 2 contributes to a renormalization 
of the nucleon's anomalous magnetic moment. In this case the renormalization is finite [48]. 

After the renormalization, the self-energy and the 'jNN vertex can be written in the 
following form: 

iZ R (p) = iX(p) - iH{M N ) - i{p - M N )E'(M N ) 

= S(p-M N ) (17) 
iT» R {p,p s ) = iY^pD-i^q^ipD+i^Fsipl) 

+i{p s -M N )^F i {p 2 s ), (18) 

where Fi^(p 2 ) = Fi^{p 2 ) — Fi^(M^) are subtracted functions. Note that functions F 3 (p 2 ) 
and F 4 (p 2 ) do not get subtracted; indeed, the Lorentz structures that correspond to these 
functions are purely off-shell — they give zero when both nucleons are on-shell (i.e., when 
both p and p s are on-shell momenta), so they do not contribute to the renormalization 
of charge or magnetic moment. However, both these functions play an important role in 
making the complete Compton scattering amplitude gauge invariant. Note also the fact that 
after the renormalization of nucleon mass, wave function, charge, and anomalous magnetic 
moment is performed, the remaining expressions for the nucleon loops (5-7) become finite. 

Now we come to the loops with A, Fig. 4. The corresponding amputated loops also give 
contributions to nucleon self-energy and to '-/NN vertex, and the corresponding expressions 
can be written in a full analogy to the case of nucleon loops. The loop (19) in Fig. 4 also 
gives a contribution to •yNN, however, this contribution is fully off-shell and momentum- 
independent: 

iT»{p,p s ) = i$fA + i(p s - M N )^B, (19) 

where A and B are constants. Nevertheless, it is important to take this contribution into 
account in order to preserve the electromagnetic gauge invariance. The renormalization of 
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nucleoli self-energy and ^NN vertex proceeds for these loops in complete analogy to the 
purely nucleon loops. 

In the case of A loops, however, we obtain in addition some higher-order divergences, 
i.e., ultraviolet divergences of &{p A ). They are to be renormalized by a corresponding ff{p A ) 
contact term. At this stage it is customary to use the MS subtraction of the higher-order 
divergences, see e.g., Ref. [68]. We have implemented the MS scheme for the higher-order 
divergences by putting the dimreg factor equal to zero (see Appendix A for more detail). 

C. Consistency with forward-scattering sum rules 

The dispersion relations enjoy a special role in nucleon Compton scattering, see Ref. [1] for 
a review. First of all, practically all up-to-date empirical values of nucleon polarizabilities are 
extracted from data with the use of a model based on dispersion relations [15, 16]. Secondly, 
in the forward kinematics, the Compton amplitude can be related to an integral over energy 
of the photoabsorption cross section, which in combination with the low-energy expansion 
yields a number of model-independent sum rules. A famous example is the Baldin sum rule: 

oo 


where the sum of polarizabilities is related to an integral of the total photoabsorption cross 
section a tot over the photon lab-frame energy v. 

In general, the forward Compton-scattering amplitude can be decomposed into two scalar 
functions of a single variable in the following way: 

Tfi(u) = e'* ■ ef(u) + iff • (e** x ?)ug(u), (21) 

where e', e are the polarization vectors of the initial and final photons, respectively, and ff 
are the Pauli spin matrices. The functions / and g are even functions of v. Using analyticity 
and the optical theorem, one can write down the following sum rules: 

oo 

}(v) = m + ^J iv .-^L_, (22a) 





where /(0) = —e 2 /M N is the Thomson amplitude and a is the doubly polarized photoab- 
sorption cross section, with the index indicating the helicity of the initial photon-nucleon 
state; cr tot = §(01/2 + 03/2)- 

These sum rules should also hold for the individual contributions of the loop graphs in 
Fig. 3. In this case the photoabsorption process is given by the Born graphs of single-pion 
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photoproduction, for which analytic expressions exist [48, 69]: 



a i/2 ' °3/2 




(7r+n) (n + n) 
1/2 + °3/2 



(7r+n) (7r+n) 



.(T°P) 
1/2 



.(T°P) 
3/2 



3 / 2 " M^z/ 3 




(23) 



where C = [e^ A M 7V /(47r^)] 2 , fi = m n /M N , s = M 2 



+ 2M N v, and 



x N = (s + M 2 N - m 2 )/2s, 
x n = (s- M 2 N + ml) /2s, 



(24) 



A = (l/2s)y/s - (M N + m v fy/s - (Mi 



JV - ^tt) 2 , 



are the fractions of nucleon and pion energy (x) and momentum (A) in the center-of-mass 
frame. 

We have verified indeed that the (renormalized) p 3 loop contributions in Fig. 3 fulfill the 
sum rules in Eq. (22) exactly for any positive v. 

It is interesting to note that the leading-order pion photoproduction amplitude, which 
enters on the right-hand side of Eq. (22), is independent of whether one uses pseudovector or 
pseudoscalar nNN coupling [69]. It essentially means that chiral symmetry of the effective 
Lagrangian plays no role at this order. The latter statement can, by means of the sum rule, 
be extended to the forward Compton amplitude at ff(p z ). On the other hand, the graphs 
(12) and (13) in Fig. 3, being the only ones beyond the pseudoscalar theory, take the sole 
role of chiral symmetry. In the forward kinematics these graphs indeed vanish but play an 
important role in the backward angles. Without them the values of a and (5 would be entirely 
different. The value of a + (5 would of course be the same, but a — f3 would (approximately) 
flip sign. Furthermore, in the chiral limit, the value of a — (3 would diverge as 1/m 2 (instead 
of l/m n as it should). We thus arrive at the conclusion that chiral symmetry of the effective 
Lagrangian plays a more prominent role in backward Compton scattering. 

D. Error due to 0(p 4 ) effects 

In the previous publication [54] we simply adopted the error estimate from Ref. [43]. 
However, in this work we compute to one order higher than in Ref. [43] and hence the error 
analysis needs to be revised accordingly. 

An error of an effective-field theory calculation is an estimate of higher-order effects 
assuming their natural size. The higher-order effects not included in our calculation begin 
at €?(p 4 ), i.e., the order at which the polarizability LECs, 5a and 5 f3, arise. The naturalness 
assumption requires these constants to be of order of unity in the units of the chiral symmetry 
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breaking scale of a GeV. To be more specific we assume the absolute value of these constants 
(in the MS scheme) is limited by 

«(err.) = P(err.) = (e 2 / An) / Mfj « 0.7 x l(T 4 fm 3 , (25) 

This number gives a natural estimate of the error on polarizability values we have obtained 
at NNLO. It is not difficult to find how this error propagates to observables, once the effect 
of the LECs on those observables is known. For example, for the unpolarized differential 
cross section, the error is given by [cf. Eq. (32)]: 

= 8 (e 2 /4vr) $ „„' [(1 + z) 2 a\ err) + (2z) 2 {3 2 err ,] 1/2 (26) 
where z = cos^iab, v [y') the laboratory energy of the incident (scattered) photon, and 



(center-of-mass frame) , 
u^(tY (lab frame) 



Let us emphasize that we do not include the errors due to the uncertainty in the values 
of parameters in Table II or due to the <ff(p 5 / A 2 ) effects which stem from graphs with two 
A propagators. Our errors are thus underestimated, however, they can directly serve as an 
indicator of sensitivity to the polarizability LECs at ff{p A ). 



IV. PROTON POLARIZABILITIES 

The chiral-loop contribution to scalar polarizabilities of the proton which arise from the 
NNLO calculation of the Compton amplitude is given in the Appendix A. In addition, we 
have the tree-level A(1232) contribution from graphs (14) in Fig. 4 and its crossed, given 
by [43]: 

a (A-excit.) = --— -|^!_ — ~ -0.1 , (28) 
v ; 47t(M 7V + Ma) 3 v ; 

P (A-excit.) = - ,J 6 % , 9 - ~ 7.1 . (29) 

Here and in what follows the numerical values are given in the units of 10 _4 fm 3 . 
The numerical composition of the full result thus looks as follows: 

a = 6.8 +(-0.1) +4.1 = 10.8, (30) 



(3 = -1.8 + 7.1 -1.3 = 4.0. (31) 

As explained earlier, a natural estimate of ^(p 4 ) contributions yields an uncertainty of at 
least ±0.7 on these values. In Fig. 5 this result, shown by the red blob, is compared with 
the empirical information, and with the A- less ^(p 4 ) HB^PT result of Beane et al. [42]. 

We can clearly see a few-sigma discrepancy of our result with the TAPS-MAMI deter- 
mination of polarizabilities [9]. On the other hand, as shown in the next section, our result 
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FIG. 5: (Color online) The scalar polarizabilities of the proton. Our result is shown by the red 
blob. The A-less HB%PT result [42] is shown by the grey blob. Experimental results are from 
Federspiel et al. [5], Zieger et al. [6], MacGibbon et al. [8], and TAPS [9]. "Sum Rule" indicates 
the Baldin sum rule constraint on a + (3 [70]. "Global average" represents the PDG summary [37]. 

agrees with TAPS data for the Compton differential cross sections. Of course we compare 
with the data at the lower energy end (below the pion threshold) where polarizabilities play 
the prominent role. The extraction of the polarizabilities in Ref. [9] has also been influenced 
by data above the A-resonance region to which we cannot compare. Clearly an extraction 
of scalar polarizabilities based on the data of 400 MeV and higher could be affected by un- 
controlled model dependencies and needs to be avoided. Excluding the higher-energy data 
from the TAPS analysis could help to resolve the apparent discrepancy between theory and 
experiment in Fig. 5. 

In Fig. 6 we show the pion mass dependence of proton polarizabilities in both B^PT 
and HB%PT. The difference between the two for the magnetic polarizability (lower panel) 
at €?(p 3 ) is stunning (compare the blue dashed and violet dotted curves). The region of 
applicability of the HB expansion is apparently limited here to essentially the chiral limit, 
0. For any finite pion mass, the B%PT and HB%PT results come out to be of a 
similar magnitude but of the opposite sign. A similar picture is observed for the irA loops 
arising at C{p A / A). In fact, we have checked that in the limit of vanishing A-nucleon mass 
splitting [A — > 0), the considered ttN and 7rA loops give (up to the spin-isospin factors) 
the same result. 

The total effect of the Delta here is the difference between the €?(p 3 ) and the ff(p A /A) 
curves in B%PT and the difference between the ff(p 3 ) and the &{e 3 ) curves in HB%PT. Note 
that here the &{e 3 ) contribution in HB^PT with A(1232) precisely corresponds to ff(p A /A) 
of B^PT, thus we do not include the iff(p A ) LECs in neither of the calculations. The actual 
C{e 3 ) calculations [44] are supplemented with the <ff{p i ) LECs, whose main role is then to 
cancel the large contribution of the A-isobar. 
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FIG. 6: (Color online) The pion-mass dependence of proton polarizabilities. Data points at the 
physical pion mass represent the PDG values. The legend for the curves is in the upper panel. 



V. RESULTS FOR CROSS SECTIONS 

In this section we present the results for the differential cross sections of proton Compton 
scattering. 



A. Unpolarized 

In Fig. 7, we consider the unpolarized differential cross section of the 7p — > jp process as 
a function of the scattering angle in center-of-mass system, at fixed incident photon energy. 
In Fig. 8, we study the same cross section, but as a function of the energy for fixed scattering 
angle in the lab frame. Our complete NNLO result is shown by red solid curves with band 
indicating the theory error estimate given in Eq. (26). The agreement between the theory 
and the experiment is quite remarkable here, especially given the fact that the theory result 
here is a prediction in the sense that has no free parameters. Despite this good agreement, 
as already noted above, there is a few-sigma discrepancy in the polarizability values between 
this theoretical prediction and the most precise empirical extraction [9]. This is apparently 



15 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — rn r~i — i — i — i — i — i — i — i — | — i — i — i — i — i — i — i — r 




FIG. 7: (Color online) Angular dependence of the 7^ — > "fp differential cross section in the center- 
of-mass system for a fixed photon-beam energies as specified for each panel. Data points are from 
SAL [7] — filled squares, and MAMI [9] — filled circles. The curves are: Klein-Nishina — dotted, 
Born graphs and WZW-anomaly — green dashed, adding the p 3 nN loop contributions of B^PT 
- blue dash-dotted. The result of adding the A contributions, i.e., the complete NNLO result, is 
shown by the red solid line with band. 



because the data at higher energies used additionally in the empirical extraction play an 
important role in the determination of (3. 

It is always interesting to study the convergence of the chiral expansion. In these figures 
the leading-order, ff(p 2 ), result is shown by dotted curve, which is nothing else than the 
Klein-Nishina cross section (i.e, Compton scattering off a classical pointlike particle with the 
charge and mass of the proton). The NLO, 0(p 3 ), result is given by the blue dash-dotted 
curve. One can see that the size of the effects varies strongly with the scattering angle. 
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FIG. 8: (Color online) Energy dependence of the "yp — > "fp differential cross section in the lab 
frame for fixed values of the scattering angle. Data points are from: Illinois [5] — open squares, 
MAMI [6] — filled triangles, SAL [8] — open diamonds, and MAMI [9] — filled circles. The legend 
for the curves is the same as in Fig. 7. 



At energies below the pion-production threshold, the NLO effects are tiny at backward 
angles but play a crucial role at forward angles. The situation is quite the opposite for 
the NNLO A-isobar contributions. Nevertheless, the convergence of this expansion seems 
to be satisfactory and in any case is much better than it would be in analogous HB^PT 
calculations. 

For completeness the result for the Born contribution, given by the Powell cross section 
together with the WZW anomaly contribution (graphs in Fig. 1), is shown here by the green 
dashed curves. Any deviation from these curves at low energies is attributed to polarizability 
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w 50 100 150 90 180 

v (MeV) 9 lab < de §) 

FIG. 9: (Color online) Linearly polarized differential cross sections as denned in the text for fixed 
values of the scattering angle (left panel) or fixed energy (right panel) . The green dashed curves — 
Born graphs and WZW-anomaly; black long-dashed — in upper panels result of adding a = 10.8, 
in lower panels result of adding /3 = 4, both using LEX. The NNLO result is shown by the red 
solid line with band. 



effects. More specifically, 
Ha ^rr (Born) 

— - — = -8(e 2 /47r) $ vv> [\{a + + zf - \{a - - zf] + (32) 

where $ is defined in Eq. (27) and z is the cosine of the lab-frame scattering angle. 

The difference between the dashed (Born) curves and the dotted (Klein-Nishina) curves 
arises mainly due to proton's anomalous magnetic moment. The difference is substantial but 
at backward angles can be seen to cancel almost entirely against the chiral loop contribution 
at 0{p 3 ), to obtain the blue curves. Thus, at ff{p i ), there is an intricate cancellation between 
the anomalous magnetic moment and the chiral loop effects. It would be interesting to see 
if this cancellation persists at higher orders. 



B. Polarized 



In Fig. 9 we show the results for Compton-scattering differential cross sections obtained 
with linearly polarized beam. The subscript x (y) indicates that the beam polarization is 
parallel (perpendicular) to the scattering plane, 9 is the scattering angle in the lab frame 
and dfl = — 2n sin 6 d6. The two particular combinations of polarized cross sections, seen 
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in the upper and the lower panels, are chosen such that each of them is sensitive to only 
one of the polarizabilities. They are therefore selected for the forthcoming measurement of 
proton polarizabilities at the HIGS facility [71]. The HIGS measurements are planned to be 
taken at 110 MeV photon lab energy, where a low-energy expansion (LEX) is assumed to 
hold. Indeed a study of the unpolarized differential cross section indicates that LEX can be 
trusted to energies of about 100 MeV [8]. Here we make a similar study for the polarized 
cross sections. 

The LEX states that at the second order the deviation of the polarized cross sections 
from the corresponding Born result (dashed lines in the figure) is simply given in terms of 
the polarizabilities: 

8(e 2 /4vr) $ vv (a cos 9 + (3) cos 9 + €?(z/ 3 ) (33a) 
8(e 2 /47r)$z/z/ / (a + (3cos9) + €?(z/ 3 ), (33b) 

where $ is defined in Eq. (27). A derivation of these expressions is given in Appendix B. 

From Eq. (33) one can indeed see that, at this order in LEX, the difference of the polarized 
cross sections, da y /dQ — dcr x /dQ, is proportional to a, while the combination cos 2 9 da y j 'dfl — 
da x /dfl is proportional to (5. One should realize, though, that this is only an approximate 
result which breaks down at sufficiently high energies. The Fig. 9 attempts to address this 
issue in a quantitative way by comparing the second-order LEX (long-dashed curves) with 
the result of NNLO B%PT (red solid curves with the error band). The LEX and B%PT 
results have exactly the same values for the polarizabilities a and f3, but the validity of 
B%PT extends over the whole considered energy range. 

We conclude that a determination at 110 MeV based on a second-order LEX can be 
reliable for a, see the upper panels. The situation is not as fortunate for the observable 
aimed at the determination of /3, see the lower panels. The LEX result begins to fail here 
at lower energies, at least in the backward angles where Compton experiments are usually 
simpler. 

VI. CONCLUSION 

We have completed a next-to-next-to-leading order (NNLO) calculation of low-energy 
Compton scattering on the proton within the ;\PT framework. More specifically, we have 
computed all the effects of order p 2 , p 3 , and p 4 /A, with A being the excitation energy of the 
A(1232) resonance. These are all predictive powers in the sense that no unknown low-energy 
constants (LECs) enter until at least one order higher [i.e., ff{p A )}. This fact together with 
the availability of precise data for Compton scattering has given us a unique opportunity to 
put %PT to a test. 

We have found that, assuming a natural size of the <^(p 4 ) LECs, the theoretical uncer- 
tainty of the NNLO calculation is comparable with the uncertainty of present empirical 
information about the cross sections of proton Compton scattering and the corresponding 
values for isoscalar polarizabilities of the proton. Within these uncertainties the NNLO 
result agrees with the cross sections data below the pion threshold but shows a three-sigma 
discrepancy in the value for the magnetic polarizability. We note that the state-of-the-art 
empirical value for the polarizabilities was extracted by using not just the low-energy data 
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but also data above the A-resonance region. The planned experiments at HIGS could be 
very helpful in sorting out this issue, since they plan to use precision low-energy data only. 
In this case, however, the reliance on the strict second-order low-energy expansion might be 
a problem, as our calculation has shown. The framework itself could provide a more 
reliable energy interpolation needed for the extraction of polarizabilities. 

In this work we have insisted on the fact that chiral power counting should be done for 
graphs, not contributions. It does not put any constraint on how many powers of pion mass 
or energy may appear in the result. It puts the constraint on the leading power only. The 
heavy-baryon expansion is therefore not mandatory for correct power-counting. What is 
important is that no powers lower than given by power-counting are present in the result. 
The manifestly covariant baryon %PT (B%PT) conforms to this requirement, because even 
if the lower-order terms appear in calculation of a given graph, they are shown to contribute 
only to a renormalization of the LECs. In our example, the low-energy theorem and chiral 
symmetry ensured that all such troublesome terms contributed only to the renormalization 
of nucleon mass, charge, and the anomalous magnetic moment. 
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APPENDIX A: CHIRAL LOOP CONTRIBUTIONS TO POLARIZABILITIES 

Hereby we give the expressions for the loop contributions of <ff{p i ) and 0(p A /A) to 
isoscalar proton polarizabilities a and f3, as well as the corresponding heavy-baryon results. 

1. Nucleon Loops 

Our results for the loop contributions at 0(p 3 ) agree with [36]: 
° = M^i - 1 + J [D^W 2x*(-3xW- te+ 5) (AD 




+ x 2 (9x 4 - 2Qx 3 + 29x 2 - 18x + 7)/T - (9x b - 33x 4 + 45x 3 - 27 x 2 + 7x - 1 )fi 

13 = 192^M N P (A2) 

2x 2 {Qx 3 -Ux 2 + 9x-l) + (9x 4 -24x 3 + 2lx 2 -6x+l)fi 2 

where \i = m n /M N , and D N {x) = ji 2 {l — x) + x 2 . 
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The corresponding heavy-baryon result is obtained from these expressions by expanding 
in jj, and keeping the leading term only: 

* {HB) = *°fff , (A3) 

76o7T _/ 777.Tr 

P iHB) = 7 J 2 9 } 2 ■ (A4) 

2. Delta Loops 

The ff(p 4 /A) loops of Fig. 4 give the following contribution to polarizabilities: 

{i 
^ + 85-3 / . dX *^ Ul - a;) (35- 104a; + 17x 2 + 112x 3 -60x 4 


+ (l05-273a;+72a; 2 +92a; 3 )5 + ( 105- 269a; +88x 2 + 72x 3 )5 2 + (35-100a;+64a; 2 )5 3 ) 
- x (35 - 69a; - 40a; 2 + 72x 3 + (35 - lOOx + 64a; 2 ) 5) li 2 ] 

-6 y"rfa;a;(l2 + 9a;-34a; 2 + 3(4-5a;)5)[S-log/JA(a;)] >, (A5) 

o ' 
( i 

e 2 h 2 A M N J 65 f dxx 2 2 







+ 6^" rfa;a;(l2 + 7a; + 10a; 2 + 3(-4 + 5a;)5)[S-logDA(a:)]|, (A6) 

where \x = m^jM^^ 5 = A/M N) and D A (x) = (1 — x)[(l + 5) 2 — x] + a;/! 2 . Furthermore, in 
these expressions we have S = 2/(4 — 6?) — 7^ + log (47r A/Mat) the divergence in d dimensions, 
with A the dimreg scale. Thus, the A loops contain an ultraviolet divergence which is to 
be renormalized by €?(p A ) LECs. We choose to define the values for these LECs in the MS 
scheme, and hence put S = 0. 

Expanding these results in small fi and 5 to leading order, we reproduce the heavy-baryon 
result for the 7rA-loop contributions [45]: 

" (HB, -8^k( 9 + '°^)' < A7 > 
p 2 h 2 9A 

o(HB) = eh A ZZ1 

P 864tt 3 / 2 Z\ g m^ { } 

APPENDIX B: LOW-ENERGY EXPANSION FOR CROSS SECTIONS 

The differential cross section is given in terms of the Compton amplitude by 

Tt = 16n(s-M 2 N y £ I TaV ' Act ' 2 ' (B1) 
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where A and a are the target and the photon's helicities. The sum is over the final helicities, 
the initial ones are fixed. To find the low-energy expansion (LEX) of this quantity at the 
second order in energy, we can ignore the spin- dependent contribution and write the Compton 
amplitude as follows: 

T\' c ',\a = (-A^s, t) S' c , ■ S c + A 2 (s, t) q ■ S' c , q' ■ S a ) 2M N 5 yx , (B2) 

where Ai are scalar amplitudes dependent on the Mandelstam variables only, q (q 1 ) is the 
initial (final) photon 4-momentum, and 

<f " = e» - ^"V , (B3) 

g'* = e' 11 - |r^<T (B4) 

with e the vectors of photon polarization, and P = p + p' the sum of the nucleon external 
momenta. Since 

Y,&? = -!T, e ff . e ; = -l, (B5) 

we have 

yj ce = _^ + ^!™_£W (B6) 

with P ■ q= i(s - M 2 N -u + M 2 N ) = M N (u + v'). We therefore obtain 

J2\ T ^'^°\ 2 = (2M N ) 2 (-A 1 ^ ll + A 2 q fl q / -^)(-A 1 g: u + A 2 q u q'-0 

AV 

x ( - nW 4- q + q - tULJL) 
[ 9 + P-q {P-q'f h 

We next take the Born contribution out of A\. 

A 1 =A 1 + e 2 /M N (B7) 

and use the LEX: 

A 1 = 47r(a + p z )vi/ + @{y z ) , (B8) 
A 2 = -AttP + ®{v) , (B9) 

where z is the cosine of the lab-frame scattering angle, a and (5 are respectively the electric 
and the magnetic polarizability. 

At the second order in v for the non-Born (NB) contribution we thus have 

E I 7 ™.! 2 = -8M N (Ane 2 )uu' \a + (3z + (q' ■ e a ) 2 (3 - j^^tf • 0> + P) 

= -8M w (47re 2 )z/z/' [(l - (<f • e CT ) 2 )a + 0z] + €?(z/ 3 ) , (BIO) 

where q' — (1, Vl — z 2 , 0, z). Substituting in Eq. (Bl) and selecting the appropriate photon 
polarization we arrive at Eq. (33). 
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